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Abstract 


In heavy ion reactions , neutron production at forward angles is observed to occur 
with a Gaussian shape that is centered near the beam energy and extends to energies 
well above that of the beam. This paper presents an abrasion-ablation model for mak- 
ing quantitative predictions of the neutron spectrum. To describe neutrons produced 
from the abrasion step of the reaction where the projectile and target overlap , we use 
the Glauber model and include effects of final-state interactions. We then use the pre- 
fragment mass distribution from abrasion with a statistical evaporation model to esti- 
mate the neutron spectrum resulting from ablation. Measurements of neutron 
production from Ne and Nb beams are compared with calculations , and good agree- 
ment is found. 


Introduction 

Neutrons produced from the nuclear interactions of 
cosmic rays with the Earth’s atmosphere, aircraft or 
spacecraft structural shielding, or the self-shielding of 
the human body are responsible for a large fraction of the 
energy deposition to passengers on high-altitude aircraft 
(ref. 1) and to astronauts in low Earth orbit. The nuclear 
reactions leading to the production of secondary neutrons 
from cosmic rays are dominated by the nucleon compo- 
nent; however, a significant fraction of neutrons are also 
produced in the interactions of the primary helium and 
heavy ion components of the cosmic rays with shielding. 
(See ref. 2.) Models of neutron production in heavy ion 
reactions are thus important inputs for the assessment of 
radiation damage from cosmic rays. 

The study of neutron production from reactions 
induced by heavy ions may also present important 
insights into the theoretical modeling of the production 
of heavy fragments in these reactions. The abrasion- 
ablation model has been used for many years to describe 
mass yields in heavy ion collisions. (See refs. 3 to 6.) 
However, few attempts have been made to calculate 
nucleon production, including momentum distributions, 
in the abrasion-ablation model (ref. 7). Heavy ion frag- 
ment mass yields and nucleon production are ultimately 
related in these reactions, theoretically through the equa- 
tions of motion or scattering amplitude. In references 8 
and 9, measurements of inclusive neutron production in 
heavy ion collisions at 390A MeV and 800A MeV sug- 
gest that neutrons from the knockout stage of abrasion 
and the evaporation stage of ablation can be separated 
from the data. Madey and coworkers (refs. 8 and 9) have 
decomposed the neutron production data at forward 
angles into three Gaussian components of increasing 
widths. The narrowest component was attributed to evap- 
oration neutrons, for which they found an effective tem- 
perature of around 2.7 MeV at 390A MeV and 3.3 MeV 
at 800A MeV. The second component of intermediate 
width was attributed to direct knockouts of neutrons. For 
this component, the width of the distribution was related 


to the internal momentum distribution of nucleons with 
Fermi momentum of about 250 MeV/c. The third and 
widest component was attributed to a high-momentum 
tail in the internal momentum distribution or to collective 
excitation effects. These results are useful because they 
suggest that some details of the different stages of the 
abrasion-ablation description can be uncovered from 
inclusive data where only the total distribution of neu- 
trons is measured. 

In this paper, we extend the abrasion-ablation model 
as formulated by Hiifner, Schaffer, and Schurmann to the 
evaluation of momentum distributions for nucleon pro- 
duction in heavy ion collisions (ref. 4). Using the 
Glauber model (ref. 10), we consider the spectrum of the 
knockouts in the overlap region of the collision. We also 
estimate the contribution of the final-state interactions of 
the knockouts for the projectile interacting with the pre- 
fragments. As in reference 4, the main approximations, 
other than the Glauber approximation, are derived from 
the treatment of the nuclear wave function, for which 
single-particle wave functions are used at all stages. In 
references 11 and 12, this formulation has been used to 
calculate proton production from l2 C and 40 Ar projec- 
tiles. Excitation energies from the abrasion stage are cal- 
culated in the geometric abrasion-ablation model (refs. 3 
and 6) and used in the classical evaporation model 
(refs. 13 and 14) to estimate the neutron spectrum that 
originates in ablation. We then have a quantitative 
approach for considering the several mechanisms 
described in references 8 and 9 and can make compari- 
sons that are consistent with the models describing heavy 
fragment yields. 

The neutrons produced in the abrasion stage have 
momentum distributions largely determined by the 
ground-state, one-body density matrix of the projectile. 
To explain the high-momentum component of the neu- 
tron spectrum, we consider recent models that account 
for correlation effects based on calculations of the 
momentum distributions rc(p). (See refs. 15 and 16.) In 
many aspects, the physics of the calculations presented 



here are similar to those contained in intranuclear cas- 
cade codes that use the Monte-Carlo method. (See 
refs. 17 and 18.) The formalism we present is useful first 
because of its simplicity because it involves only a few 
numerical integrations and second because of its ability 
to test nuclear structure inputs, such as the one-body 
density matrix. In the remainder of this paper, we first 
introduce the Glauber amplitude and recast the abrasion 
model in terms of the momentum distributions of the 
knockouts. The final-state interactions are then studied 
by correcting transition densities for rescattering effects. 
The energy spectrum of neutrons decaying from pre- 
fragments are then considered by using the classical 
evaporation model. Finally, we make comparisons 
with experiments and discuss the results of model 
calculations. 


Symbols 


A 

B 

b 

c 

E 

F 

F* 

F i 

f 

Sn 

Im 

y (1) 

m 

K 

k 

M 

m N 

N l 

n 

n(p) 

«(x) 


P 

*( b, b') 
Q 


mass number 

two-body slope parameter, fm 2 

impact parameter, fm 

speed of light, m/sec 

energy, MeV 

fragment 

prefragment 

probability of emission of ion l 

scattering operator, fm 

statistical weight 

imaginary part of function 

cylindrical Bessel function of first kind of 
order m 

projectile target relative wave number, fm -1 

wave number of emitted neutrons, fm~ 3 

mass, MeV/c 

nucleon mass, MeV/c 2 

single-collision term 

number of abraded nucleons 

momentum distribution 

Fourier transform of nucleon momentum 
distributions 

projectile 

function describing projectile spectators 
defined in equation (13) 


q momentum transfer, frrT 3 

r internal nuclear coordinate, fm 3 

S n separation energy, MeV 

s transverse part of r 

T target 

level density of residual nucleus 
X final target state 

z component of r 

a ratio of real to imaginary parts of f NN 

P relative projectile target velocity 

T profile function 

8 Dirac delta function 

£ target constituent 

u neutron mass 

* n 

% defined in equation (22) 

A n defined in equation (20) 

p(r) one-body density, fm -3 

p (r, r) one-body density matrix, fnT 3 
a cross section, mb 

<]>(r) single-particle wave function 

X Eikonal phase 

¥ complete nuclear wave function 

ft Eikonal inelastic collision term 


Subscripts: 

abl ablation 

abr abrasion 

CN formation cross section 

/ final state 

i initial state 

j abraded nucleons (projectile constituents) 

NN nucleon-nucleon (two body) 

P projectile 

T target 

X final target state 
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Abrasion Theory 

In the Glauber model the scattering operator for nucleus-nucleus collisions, as given in references 10 and 1 1, is 

f=l 2 1 p 2fcexp(iq b)r(b) 0) 

where K is the relative wave number of the projectile target, b is the impact parameter, and q is the momentum transfer. 
The profile function representing the multiple-scattering series at high energies is 

r(b) = i-n [ 1 " r cy(b- 8 «-S/)] ( 2 ) 

u 

where £ andj label the target and projectile constituents, respectively. In equation (2), is the two-body profile func- 
tion with the internal coordinate having components r = (s, z). 

The scattering amplitude of equation (1) is related to the production cross section for a projectile nucleon from the 
abrasion process by 


do 
d k 


= Y— !— f dE .d 2 qd 2 bd 2 b' exp [/q • (b-b')] 8(£.-£ f ) 
^(2n) 2 J F 1 1 


m 

j = 2 


n r dk. 1 , + ,| 

— 1 - 1 (TPlr (b ) |XF*k.)(k.F X|T (b) \PT) 


L(2tc) 3 J 


r ' j 


(3) 


where k y denotes the wave numbers of the abraded nucleons, F* denotes the prefragment (with A f . = A p -n), and 
we have inserted initial and final states in equation (3). 

The state dependence of the final target energy prevents closure on these states from being automatic in 
equation (3); however, when energy conservation is not considered, it is made outright (ref. 4). Next, we consider the 
change in energy of the target from the collision as follows: 


Ej-E^ = Ej- JiPj-- q) 2 + Afjf = E \ 1 - 1 + 




Ej 


(4) 


where Mj and M x are the mass of the target in the initial and final states, respectively. We expect that performing clo- 
sure on the IX) will be valid for sufficiently large values of Ej.. After closure on the final target states, we find 


S = ^2 j dE p .d 2 qd^bd 2 b' exp [iq ■ (b -b') ] C^b, b', k, q, E p ?j 


where we define 


C J1 (b,b',k,q,E r ) = <71 jin '~~3 5(E j .-£ / )<F|r^b')|F*k.)<k.F*|r(b)|F)||7) 


j = 2 


To consider the energy-conserving 6 function in equation (6), we introduce the Fourier transform pair 

o n (t) = jdE exp (-iEt) o n (£) 


and 


O n (E) = exp (iEt) o n (t) 


(5) 


( 6 ) 


(7) 

( 8 ) 
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In the projectile rest frame, we have 


<» 

/ = I N 

where S n is the separation energy and T is the recoil energy of the prefragment, including any excitation energy of 
the prefragments. We next go into temporal space to consider the dk . integrals in equation (6) and thus rewrite 
equation (9) as 


From equations (6), (7), and (10), we find 


dk. 


=<*i m 

j = 2 


L (2tc) ■ 


_ k? 
E.-E f = £-Y-i. 
' f A* 2m w 
j = 2 ■ « 


k?r 


( 10 ) 


exp X 2 ^" | (P|r (b<) I F’k.Xk.nr (b) |P) |7) 


j-Ti 2m N 


( 11 ) 


To simplify equation (11), we first factor the profile function into projectile participant and spectator terms as 

Ap n 

r(b) = 1- [J JJfi.(b-s.) (12) 


l = n+\ 


7=1 


where 


Oj- no-y 


(13) 


In the abrasion model, the orbits of the prefragments are assumed to be nearly the same as those of the projectile. 
This assumption is consistent with the use of the impulse and frozen nucleus approximations at high energies. A com- 
pletely factored form in the participant and spectator coordinates is assumed for the projectile wave function 


IP) = 10140 


(14) 


where | P) and |<|> n ) are the wave functions of the core (spectators) and of the knockouts (participants), respectively. The 
antisymmeterization is ignored in equation (14), which should be accurate if the mass of Pis much larger than that of the 
knockouts. Antisymmeterization in the subsystems of |F) and |<t>„) may still be included. By using plane-wave states for 
the |kp> and substituting equations (12) and (14) into equation (1 1), we find that 

r f A ) 

o„(') = (7| | F (PinC/lb'-sjllP’XF'lTTe^b-s^lpJrfr^r' exp(ik x,)(2^(b'-s')2(b-s 1 ) 

[ v"/ ; 

n 

*11 

7=2 

where xj = r j - r j . Using the coherent approximation for the target wave function in the intermediate states and the 
independent panicle model for the fragment wave function leads to the following equation: 


dk 


f i 

J(2 exp 


r ~k 2 t 


♦:cr; r„) M7) 


(15) 


G n {t) 


(A ' 
A p 

K n J 


dN, 


? A r- n (b, b ) A n _ , (b, b , t)-^ 


(16) 
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where the function !P(b, b ) describes the projectile spectators as given by 

A“ n (b,b') = mUQ^ - s ;)nr iFte (b - »,> ro 


( 17 ) 


i 


i 


We next perform closure on the prefragments states in equation (17) because we do not consider coincidences with indi- 
vidual states. After closure we find that 


'P Ap (b,b') = <rF|f[G^b'-s;]e z (b-s ; )|F7) 


(18) 


In equation (16) we have defined 

~d£ ~ ~(^TCp^ dTdT exp ^ ,k x ^P^ r ‘ r ^/ + ( b ~ s )Qi(b-s) 


(19) 


where p(r, r ) is the one-body density matrix of the projectile given by p(r, r) = <j> + (r ) <)> (r) . Next from 
equation (15), after evaluation of the integrals over k . for j > 2, we find that 


A J| _ 1 (b, b', t) = 

7 = 2 


dr. dr. 


m 


N 


J J\2nit 


\3/2 

j exp| 


f r\ 

m Nj 

lit 


p(r 7 ,r;)e/( b '- s ;)e ; .(b-s.) 


I D 


( 20 ) 


In energy space, 

Vl( b ’ b '’V k ) = <71 J n [ Jr A :p ( r / r j] Q j{ b " S 'j) Q j (b “ s 7 } ] 


7 = 2 


m 


N 


2 \2n 


, N ,, n/ ,M3(»-l)/ 2 l-l 

1 A3 (n- l)/2 S„_ i 


X l 3 1” 
n - 1 


- l)/2] - 1 J {Wn-\)/2\-\ ^n-i X n-0^ 


( 21 ) 


where J ( 1 ) is the cylindrical Bessel function of the first kind of order m and where 


E , = 2m 
- i 


N 


k 2 

T + K 

F 2m 


N 


~ S n~ ( E T~ E x) 


( 22 ) 


and 




(23) 


For n - 1 , we have Aq = 5 ( T F ,). If we assume forward-peaked density matrices (about), a small argument expansion 
of the Bessel functions can be developed (ref. 11) that results in 




(«-0 


■ 5„_, 1 




where, for example, C x = 1, C 2 = C 3 = y^:, and C 4 = 

The nucleon momentum distribution from abrasion then takes the following form: 


Jn- 1 


(24) 


(ML - (b-b')J A-"(b,b)^J</7- F .A„.,(b,b;7- r ,k) (25) 


( 271 ) 
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If 


q 1 + (26) 

we may approximate equation (25) as 

(sL ■ , <b, v k> a7) 

The result of equation (27) suggests that for A p » 1, the momentum dependence of higher production terms (n > 1) 
should be similar to that of the leading-order terms. (See appendix.) This result supports the success of the hard- 
scattering model of Hatch and Koonin (ref. 19), which only uses the single-scattering mechanism to predict the shape of 
the inclusive proton distribution in heavy ion collisions. The model developed here differs from the hard-scattering 
model by our use of the target closure approximation. In this approximation, the effects of smearing of the secondary 
momentum from the target knockouts are not considered. Instead, these effects are replaced by averages represented by 
the target density and by using only on-shell, two-body amplitudes. Also, the Glauber model has a much fuller multiple- 
scattering structure than does the hard-scattering model. 

To include the effects of final-state interactions of the nucleon knockouts, we use the Eikonal model described in 
reference 1 1. In this model, the plane waves are replaced by the distorted wave for the nucleon-projectile recoil interac- 
tion evaluated at the relative energy between the knockout and the recoil. Modifying equation (19) as in reference 1 1 
gives 


aiv i i r , 

~dk = J2k)3J ] ^ ex P( /k ' x )P( r ’ r ) exp[-2/mx (_) (y)]Qj (b' -s') Q^b-s) (28) 

where % ^ is the outgoing Eikonal phase. Equation (28) ignores off-shell effects but includes the energy dependence of 
the final-state interaction and assumes a medium modified interaction, as described in reference 11. 

Optical Limit for Profile Functions 

For ApAj» 1, the optical limit of the profile functions occurring in the previous equations may be used (refs. 10 
and 20). From reference 20, we find in the optical limit that 

<P Af, (b, b') = exp { i [x (b) - x + (b') ] } (29) 

where the Eikonal phase is 


A pAj » 

x(b) = J2nk^r q exp ( iqb ) F p^) F r^)/yv/v^) (30) 

with F (q) denoting the one-body form factor and f NN denoting the two-body scattering amplitude, which we represent 
by 


?NN ~ 




4?t 


AW f lf> 

— exp^yVJ 


(31) 


where is the two-body total cross section, is the relative momentum in two-body center of mass frame, B the 
two-body slope parameter, and a the ratio of the real to imaginary parts of f NN ( q = 0). For the inelastic terms, we write 
in the optical limit (ref. 20) that 

b ' - SjjQj (b - Sj) = exp [£2^ b' - s ', b - sj] - 1 (32) 

with 

Q(V- S :,b-s.) = 2 XPVV exp [iq ■ (b-s. + s a )j exp[-iq • (V - s’ + sj]/^ (q)f NN (q)(33) 
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Model for Nuclear Density Matrix 

We next describe a local density approximation for the one-body density matrix. For a projectile nucleus, the one- 
body density matrix is defined in terms of the complete nuclear wave function 'F (ref. 15) as follows: 

p (r, r) = jdr 2 dt^ ... ^'F+jV, r 2 , ... , r A J'FjV r 2 , ... , r A J (34) 


Evaluation of equation (34) requires knowledge of the complete nuclear wave function; however, in practice a model is 
introduced. In the Fermi gas model, the density matrix is represented as 


P(r, r) 


3y 1 (fc F lr-rl) 


k F \r-r | 


(35) 


where kp is the Fermi momentum and p g is the density of nuclear matter. The Fermi gas model is known to provide a 
poor representation of the density matrix. However, its form suggests the use of a local density model, for which the 
density matrix is assumed to factor as 


p(r, r ) » p (y)n(x) 


(36) 


with x = r - 1-’ and y = - (r + r ) and where the one-body density is given by the diagonal part of the density matrix 


P (r) = p (r, r ) (r = r) 

and n (x) is the Fourier transform of the nucleon momentum distribution 

n(\) = jdp exp (ip • x)n(p) 

where n (p) is defined by 

n (p) = jdrdr' exp (-ip • x) p (r, r ) 

with normalization 

L (p )_f*P_ = i 

J (2n ) 3 


(37) 


(38) 


(39) 


(40) 


The one-body density is reasonably well known from elastic electron scattering. The nucleon momentum distribution at 
small to modest values of p is known from inclusive inelastic electron scattering. For large p values, backward produc- 
tion of protons suggests that large enhancements to the nucleon momentum distribution result from correlation effects 
when compared with predictions of independent-particle models. In addition, this enhancement is largely independent of 
nuclear mass. 

Haneishi and Fujita (ref. 21) have introduced the following momentum distribution: 


*(P) 



C. exp 


2P, 2 


i J 


(41) 


where no is a normalization constant. The last term in equation (41) is expected to directly reflect the nuclear correla- 
tions. The value of C 3 = 0.003, used in reference 21, was based on estimates from backwards proton production. This 
value was found to be too large in references 1 1 and 12; thus, the value Cj = 0.0008 is used herein. In equation (41), 
p j is related to the Fermi momentum by p j = j2/5kp . Values for kp are listed in table 1, and the other parameters of 
equation (41) are listed in table 2. Figure 1 compares the model of equation (39) with a Fermi gas model that includes 
only the first term in equation (41) and with the correlation model of reference 16. 


7 



Table 1. Experimental Determination of Fermi 
Momentum for Several Nuclei 


Nucleus 

kp , MeV/c 

12 c 

184 

20 Ne 

230 

40 Ar 

240 

93 Nb 

255 

197 Au 

265 


Table 2. Parameters for Momentum Distribution Model 
of Equation (41) 


i 

Pi , MeV/c 

c 

1 

jl/5k F 

1 

2 

JUlk F 

.03 

3 

500 

0.0008 



Figure 1 . Internal nucleon momentum distribution for Ne versus 
momentum for several models. 


Neutron Spectrum From Ablation 

After the projectile-target collision, a distribution of 
prefragment nuclei in states of excitation remain. The 
excitation spectrum is treated as an average state in the 
abrasion-ablation model, with a single or small number 
of average excitation energies used to describe the pre- 
fragment and the strength of the state determined by the 
total abrasion cross section. A microscopic formulation 


of nucleon and cluster abrasion describes transitions to 
individual levels of the prefragment nuclei (ref. 22). Here 
we follow the approach described in references 3 and 6 
for estimating average excitation energies and calculate 
the neutron spectrum from ablation with the Weisskopf- 
Ewing statistical decay model. We define P n (j, k) as 
the probability that a prefragment labeled j with mass 
number Aj y charge number Zy, and excitation energy 
Ej, emits a neutron of momentum k (refs. 13 and 14). 
The momentum distribution for neutron production from 
ablation in the projectile rest frame is written as 

= <42) 
j 

The total abrasion cross section in equation (42) is evalu- 
ated from the abrasion momentum distribution in equa- 
tion (27), as discussed in the appendix. 

In the statistical model, the prefragment (compound 
nucleus) is assumed to be infinitely heavy, and the emis- 
sion spectrum is assumed to be isotropic in the rest 
frame. The probability function, as given in references 13 
and 14, is 


P nU,E n ) 


2H, 


nSn E n a CN W o 


X £ , 



(43) 


where \i n is the neutron mass, g n is the statistical 
weight, <5 cn is the formation cross section by the 
inverse process, and w o is the level density of the resid- 
ual nucleus. In equation (43) we have 


F i = 



O', E) dE 


(44) 


We consider competition between the emission of the 
light particles n, p y d y r, A, and a. The model parameters 
are taken from Dostrovsky (ref. 23). If sufficient excita- 
tion energy is available, then several neutrons will be 
emitted. We then evaluate the cumulative spectrum as 

P n^E n ) =P n U,E n ) 

+ y f 6 ' p, U, E,) P n (k, EJ dE t + ... (45) 

i J 0 


In this work, terms through the third order in the series of 
equation (45) are considered. The neutron momentum 
distribution in the lab frame is found by multiplying 
equation (42) by the neutron energy to form a Lorentz 
invariant and then performing the transformation to the 
laboratory system. 
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Results for Forward Neutron Production 


The production of secondary neutrons from 390A 
MeV and 800 A MeV 20 Ne and 93 Nb beams has been 
measured in the forward direction with targets of NaF 
and Pb (refs. 8 and 9). In figures 2 to 4, we compare these 
data with calculations based on the previous formalism. 
For calculations, the NaF target is represented by 20 Ne. 
The measurements are for the inclusive production of 
neutrons, and the calculations shown consider the neu- 
tron production from the target nuclei in the knockout 
stage of abrasion. For the Nb projectile, some under- 
estimation of the evaporation peak centered at the beam 
energy (fig. 4) occurred and may be due to the inclusion 
of only up to third-order terms in the evaporation chain. 
The underestimation of the evaporation peak for the 
93 Nb beam may not be of concern for cosmic ray studies 
for which Fe is usually the heaviest ion of interest. 

The high-energy component of neutrons extends 
well above the velocity corresponding to that of the pro- 
jectile beam. As discussed previously, in predicting 
inclusive proton yields (ref. 12), model calculations are 
only successful in predicting the high-energy component 
of the spectrum when the internal nucleon momentum 



E, MeV 

Figure 2. Calculations versus experimental data (from refs. 8 
and 9) for energy spectrum of secondary neutrons at 0° in labo- 
ratory for Ne beams at 390A MeV colliding with NaF targets. 


Neutrons from ablation of F* 
Neutrons from abrasion stage 
Total 



Figure 3. Calculations versus experimental data (from refs. 8 and 
9) for energy spectrum of secondary neutrons at 0° in laboratory 
for Ne beams at 8CKM MeV colliding with NaF and Pb targets. 



E, MeV 

Figure 4. Calculations versus experimental data (from refs. 8 and 
9) for energy spectrum of secondary neutrons at 0° in laboratory 
for Nb beams at 80(M MeV colliding with Nb targets. 
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spectrum contains correlation effects. Here the Fermi 
gas model underestimates the data by several orders of 
magnitude. The dependence on the target mass of the 
high-energy component (fig. 3) indicates that multiple 
scattering also affects this component. The model calcu- 
lations do not account for higher order cascade effects of 
projectile knockouts interacting with the projectile pre- 
fragments (other than the overall absorption). The calcu- 
lations also do not account for the production of nucleons 
through the decay of nucleon isobars that are produced in 
the reaction. Both of these contributions are important 
for energies intermediate between that of the initial beam 
and target and grow in importance with the angle of pro- 
duction. Calculations of these processes will be dis- 
cussed in future work. 


Concluding Remarks 

The production of secondary neutrons in heavy ion 
collisions has been formulated with the two-stage 
abrasion-ablation model. Good agreement with experi- 
mental data for forward neutron production is found. 
Calculations show that neutrons produced in both stages 
of the reaction are important. The secondary neutrons 
produced in nuclear abrasion extend to velocities well 
above that of the beam velocity in these reactions. This 
extension is a result of correlation effects that occur in 
the internal nucleon momentum distribution and through 
multiple scattering. The results of these calculations will 
be used to develop databases for cosmic ray shielding 
studies. 
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Appendix 

Kinematical Phase Space and Inclusive 
Momentum Distribution 


n 

dc = (2jc) 4 [dp r ] X n [dq y ] I T fi \ 2 

X r = 0 n = \j _ i 

where the phase space factor is defined as 


In this appendix, we discuss the relationship of the 
multiple-production terms of nucleons abraded from the 
projectile to the kinematical phase space. We also con- 
sider an approximation with the Glauber multiple- 
scattering series, in which energy conservation is ignored 
entirely. 


1 1 

% _ A / _ I * ~ (^4) 

P P COS 0 ^-^ L j 

v j J 

The momentum distribution of nucleons from abrasion is 
then given as 


The scattering amplitude for the heavy ion collision 
is related to the cross section by the phase space of each 
particle that appears in the final state. We consider inclu- 
sive reactions in which a nucleon originating in the pro- 
jectile is measured. For simplicity, the final target state is 
not considered, and we use closure on these states with a 
single momentum vector denoted by p-^ to represent 
these states. The cross section is then determined by 

do = i2n^Y d d Vx dp F , £ f\ [d Pj ] 

H X n = lj = 1 

x8(E r E f )8(p r p f )\T f .\ 2 (Al) 

where P is the relative projectile target velocity, F * rep- 
resents the prefragments, n is the number of nucleons 
knocked out of the projectile in the overlap region with 
the target, and i and /denote the initial and final states, 
respectively. The prefragment decays through particle 
emission when sufficient energy is available. To include 
the phase space of decay products of F*, we write 


= w'lfaT'L in*, nwn 2 <■**> 

F X n=l j =2 r =0 

Equation (A4) corresponds closely to equation (3) if we 
make the following replacement: 

K = jdq L b (£. - E f ) jdE F . 5 (£. - E f ) (A6) 

If energy conservation is ignored entirely, then we would 
have K = 1. When the participant-spectator arrange- 
ment of the Glauber series previously discussed is used, 
the inclusive momentum distribution becomes 



where 


dN 

n 

dk 


dN, 

-jrill-S'Cb)]"- 1 


and note that 


(A7) 


(A8) 


dp F . = dp F Y\dp r (A2) 

r = 0 

where r denotes the ions (if any) emitted in ablation. 
To consider nucleon production from ablation, we 
would study the p r . We use the momentum-conserving 
8 function in equation (Al) to eliminate p F from 
equation (Al). 

We next consider using energy conservation in equa- 
tion (Al). Working in the projectile rest frame, we trans- 
form dp x to dq, where q = P^-P^ is the total 
momentum transfer in the collision and use the energy- 
conserving 8 function in equation (Al) to eliminate dq L 
where q L is the longitudinal momentum transfer. We 
then find 


dN, 

J dk 


dk = 1 -2>(b) 


(A9) 


The result of equation (A8) indicates that the inclusive 
momentum distribution from abrasion essentially follows 
the shape of the leading-order term in the approximations 
discussed because since the absorptive factors in equa- 
tion (A7) change slowly with increasing n for A p » 1. 
The total cross section for abrading n nucleons follows 
from equations (A7) through (A9) as 


(A t 


C n = 


V n 


jd 2 b(P‘ 


Ap ~ n (b) [1 -2>(b)] n 


(A10) 


which is in agreement with references 4 or 5. 
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